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Abstract 

Phonon softening phenomena of the F point optical modes including the lon- 
gitudinal optical mode, transverse optical mode and radial breathing mode in 
"metallic" single wall carbon nanotubes are reviewed from a theoretical point 
of view. The effect of the curvature-induced mini-energy gap on the phonon 
softening which depends on the Fermi energy and chirality of the nanotube 
is the main subject of this article. We adopt an effective-mass model with 
a deformation-induced gauge field which provides us with a unified way to 
discuss the curvature effect and the electron-phonon interaction. 

Key words: carbon nanotube, graphene, Raman G band, phonon 
self-energy, curvature effect, energy gap, Fermi energy 



Preprint submitted to Physica E 



May 3, 2010 



1. Introduction 



The lattice structure of a single wall carbon nanotube (SW NT) can be 
specifie d uniquely by the chirahty defined by two integers (n, m) (Saito et al. 



1992alJbr)]. and the chi r ality can be deter r nined by Raman spectroscopy 



Jorio et all (12001L l2003f ): iDresselhaus et all (12005f l]. A simple tight-binding 



model shows that a SWNT is primarily metallic if n — m is a multiple of 3 
or semiconducting otherwise. A "metallic" SWNT can have a mini-energy 
band gap due to the curvature of a SWNT which gives rise to a hybridiza- 
tion between the a and tt orbitals. The presence of an energy band gap 
in a metallic SWNT has attracted much a t tention since th e early stages 
of nanotube research Hamada et al.l ( 119921 ): iMintmire et al.l ( 119921 )] . The 
present paper deals with the effect of curvature on the Raman spectra for 
two in-plan e F point lon g itudinal and transv e rse op tical phonon (LO and 



TO) modes [Farhat et al.l (12007): ISasaki et al.l (l2008bl)] and the out-of-pl ane 
radial breathing mode (RBM) [Farhat et al.l (l2009l ): ISasaki et al.l (l2008ah ]. 



In the Raman spectra of a SWNT, the LO and TO phonon modes at the 
F point in the two-dimensional Brillouin zone (2D BZ), which are degener- 
ate in graphite and graphene, s p lit int o two peaks, denoted by and G~ 
peaks, respectively, [Jorio et al.l ( 120031 ): ISaito et al.l ( 119981 . |2003|)] because of 
the curvature effect. The splitting of the two peaks for SWNTs is inversely 
proportional to the square of the diameter, dt, of SWNTs due to the cur- 
vature effect, in which does not c hange with changing dt, but the G~ 



frequency decreases with decreasing dt [Jorio et al.l ( 12002| )]. In particular, for 
metallic SWNTs, the G~ peaks appear at a lower frequenc y than the G~ 



peaks for semiconducting SWNTs with a similar diameter [Pimenta et al. 
( 1l998l )]. The spectra of G~ for metallic SWNTs show a much larger spectral 
width than that for semiconducting SWNTs. 

It has been widely accepted that the frequency shift of the (j-ba nd in 



metal li c SWNTs is produced by t h e electron-phonon (el-ph) i n teraction Piscanec et al 



2004) 
( l2006h 



Lazzeri and Mauril (20061) :llshikawa and Andol (120061 ) : lPopov and Lambin 



Caudal et al.l (|2007| ): iDas et al.l ( 120071 )]. An optical phonon changes 



into an electron-hole pair as an intermediate state by the el-ph interaction. 
This process is responsible for the phonon self-energy. The phonon self- 
energy is sensitive to the structure of the Fermi surface KohnI (119591 )] or 
the Fermi energy, E^- In the case of graphite intercalation compounds in 
which the charge transfer of an electron from a dopant to the graphite layer 



can be controlled by the doping atom and its concentration, lEklund et al. 
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(119771 ) observed a shift of the G-band frequency with an increase of the spec- 
tral width. In this case the frequency shifted spectra show that not only 
the LO mode but also the TO mode is shifted in the same fashion by a 
dopant. For a graphene mono-layer, Lazzer i et al. calculated t h e .Ef depen- 
dence of the shift of the G-band frequency Lazzeri and Mauri! (20061 ) ]. Th e 



LO mode softening in m etallic SWNTs was shown by 



Dubav et all (l2002h : 



Dubav and Kressd (20031 ) on the basis of deii sity functional theory. Recently 



Nguyen et al.l (120071 ) and lFarhat et al.l (120071 ) observed the phonon softening 



effect of SWNTs experimentally as a function of Ep by electro-chemical dop- 
ing, and their result s clear l y sho w that the LO phonon modes become soft 
as a function of Ep. lAndd (120081 ) discussed the phonon softening for metal- 
lic SWNTs as a function of the Ep position, in which the phonon softening 
occurs for the LO phonon mode. In this paper, we consider the effect of 
a curvature- induced mini-energy gap on the frequency of the LO, TO, and 
RBM in "metallic" SWNTs. 

The organization of the paper is as follows. In Sec. [2]we show that the cur- 
vature of a SWNT gives rise to a hybridization between the a and vr orbitals. 
Then we show our calculated result for the curvature-induced mini-energy 
gap appearing in "metallic" SWNTs. The current status of the scanning 
tunneling spectroscopy experimental results is briefly mentioned, confirming 
the curvature- induced mini-energy gap. In Sec. [3] we formulate the phonon 
self-energy which is given by the electron-hole pair creation process. The 
Fermi energy dependence of the self-energy is shown for graphene with or 
without an energy gap, as a simple example. In Sec. H] we provide a theo- 
retical framework for including a lattice deformation into an effective-mass 
Hamiltonian. A lattice deformation is represented by a deformation-induced 
gauge field which is shown to be a useful idea to discuss both the appear- 
ance of the curvature- induced mini-energy gap and also the el-ph interaction. 
Sec. [5] is a main section in this article in which we discuss the effect of cur- 
vature on the phonon self energy. In Sec. [6] we discuss and summarize our 
results. 



2. Curvature Effect 

Let us start to discuss the effect of the curvature of a SWNT on the 
hybridization between the a and vr orbitals (Sec. 12. ip . and we then show 
the calculated result of the curvature-induced mini-energy gap appearing in 
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"metallic" SWNTs (Sec. l2.2( l. The phonon softening phenomena are sensitive 
to this mini-energy gap. 



2.1. Curvature- Induced Hybridization 

At each carbon atom located at r on the surface of a SWNT, we define 
the atom-specific {x, y, 2;)-coordinate axes and the unit vector for each axis 
by ej(r) {i G {x, y, z}), where ez{r) is taken as the unit normal vector to the 
cylindrical surface, and ea,(r) and ey{r) are unit vectors in the tangent plane 
[see Fig. dJ^a)]. Here, e2;(r) is taken to be parallel to the axis of a SWNT. In 
the case of a flat graphene sheet, we can set the common axis vector for all 
carbon atoms and thus a unit vector ej at ri can be taken orthogonal to the 
other Bj at T2 so that ej(ri)-ej(r2) = 6ij. For SWNTs, however the orthogonal 
conditions are not satisfied because of the atom specific coordinate, that is, 
e^(ri) ■ e^(r2) ^ 1, e^(ri) ■ ey{r2) ^ 0, etc. 

To see the curvature effect more clearly, it is useful to project ei(ri) and 
ej{r2) into 

e,(ri) = e,^(ri) + e^n), e,(r2) = ef{r2) + e;(r2), (1) 

where || (_L) denotes the vector which is parallel (perpendicular) to the dis- 
placement vector r2 — ri [see Fig. [11(b)]. Let \pi{r)) {i G {x^y^z}) be the 
2pj-orbital of a carbon atom located at r. Then, the transfer integral from 
\piijci)) to \pj{r2)) may be written as 

{Pj{v2)\'H\p^{v{)) = npp^ef ir2) ■ e,^(ri) + np,Jlir2) ■ ef (n), (2) 

where Uppn and Tippa are the transfer integrals for vr and a bonds, respec- 
tively. According to a fi rst-principles calculat ion with the local density ap- 



proximation obtained by iPorezag et al.l (jl995l ). Hppn ~ —3 eV and Tip ~ 



^ppcr 



eV for nearest-neighbor carbon sites. Using Eq. ([T]), we eliminate ej-(ri) and 
ej-(r2) from Eq. ([2]), and get 

{pj{r2)\'H\pi{ri)) = Hpp^ej (r2) ■ ei(ri) + {Hpp^ - Hpp^) e5(r2) ■ e|(ri), (3) 

where we have used e^(ri) ■ e^j{r2) = and e|'(ri) ■ e^(r2) = 0. The last 
term of Eq. ([3]) corresponds to the curvature effect of a SWNT. Note that the 
coefficient of the last term includes Tippa showing that the a bond i s part ially 
incorporated by the curvature- induced hybridization [See Ando ( 2000l ) for 
more details]. 
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Figure 1: (a) The curvature-induced hybridization between two Pz orbitals of carbon 
atoms at ri and Y2 is ihustrated. ei(r) (i S {x,y,z}) denotes the basis of the {x,y,z)~ 
coordinate system whose origin is located at a carbon atom r. (b) e| (r2) and e| (ri) induce 
the hybridization, including a bonding. 



In the case of a flat graphene, we liave e2(ri) = and ez(r2) = 0. Then, 
the last term of Eq. ([3]) disappears and the theoretical model taking only the 
2pz orbital (or vr-orbital) into account becomes a good approximation. The 
curvature of a SWNT results in e!l(r2) ■ el(ri) 7^ 0, e|(r2) ■ eil(ri) 7^ 0, etc., 
and the last term of Eq. ([3]) is non- vanishing and consequently the curvature- 
induced hybridization occurs. The curvature- induced hybridization is rele- 
vant to the following two physical properties. First, the hybridization can 
open a mini-gap (up to ~ lOOmeV) near the Fermi energy in metallic SWNTs. 
Second, the curvature-induced gap depends on the SWNT (ra, m) chirality. 
For example, the gap is zero for armchair SWNTs, while it is about 70 
meV for a (12, 0) metallic zigzag SWNT. The chirality dependent curvature- 
induced energy gap will be analytically given in the next subsection. 



2.2. Curvature-Induced Mini-Energy Gap 

In Fig.|2]^a) we plot the calculated curvature- induced energy gap, -Egap, for 
each (n, m) for metallic SWNTs as a function of the chiral angle 9{°) and tube 
diameter (ij(nm). We performed the energy band struct ure calculation in 



an ext ended tight-binding (ETB) framework developed by ISamsonidze et al 
( 120041 ) to obtain -Egap- In the ETB framework, 2s and 2p orbitals, and 
their transfer and overlap integrals up to fourth nearest neighbor atoms are 
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taken into account [see iPopovl (12004 ): ISamsonidze et al.l ( l2004j ) for more de- 
tails]. We have adopted the values of the transfer and overlap integrals 
as a function of the car bon-carbon inter-atomic distance that were derived 
bv lPorezag elaP f ll995h .[^ 

Figure |2](a) shows that, for a fixed diameter of a metallic SWNT dt-, a 
zigzag SWNT {9 = 0°) has the largest value of -Egap and an armchair SWNT 
{9 = 30°) has no energy gap. The calculated results are well reproduced by 



ap 



cos 3^, 



(4) 



with c = 60(eV-nm^) [Sasaki et al.l (l2008al )]. The chirality and d i amete r 
dependence o££'gap is consistent with the results by iKane and Meld (119971 ) , 
an dlAndol (l2000h. The v alue of c is about two times larger than the result 
by iKane and Meld ( 119971 ) . This difference may come from the inclusion of 

n 



vppa in our calculation. As we will explain in detail in Sec. 14. 2[ the curvature 
moves the Dirac point in fc-space away from the hexagonal corner of the first 
BZ. As a result, the curvature can cause the quantized transverse electron 
wave vector (the cutting line) to miss the Dirac point and make a gap [see 
the inset in Fig. EJ^a)]. 

When we discuss the phonon softening of the RBM, the relationship be- 
tween the mini-energy gap and the RBM phonon energy will be important. 
In Fig. El^b), we plot the energy of the RBM, 



RBM 



Cl 



+ C2, 



(5) 



as a solid curve for comparison. Here fyjJ-Bs,M is a monotonic function of 
the tube diameter ((if[nm]) and is modeled as being linear in the inverse 
diameter, with an offset C2 which is known as the effect of the substrate. We 



^In the ETB program, we numerically solve the energy eigenequation, Hj^E') = -E'l^'), 
in the basis of |s(r)) and \pi{v)) for two carbon atoms (A and B). The basis orbitals for 
the A-atom are non-orthogonal to those for the B-atom due to the curvature effect, and 
the Hamiltonian and overlap matrices are 8x8 matrices. We assumed the on-site energies 
E{2p) = -4.882[eV] and E( 2s) = -13.573feV]. £;(2s) - E{2p) « -8.7[eV] is close to the 
value (-8.868[eV]) shown in lSaito et all (|l992bl) . 

^ Although the energy gap at the Fermi level has little to do with the overl ap integral, we 



shall n ote that the overlap integrals S*^*^ and S*^^ are switched in Table 1 of lPorezag et al 
(|l995h . 
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0.5 1 1.5 2 2.5 



dt{nm) 



Figure 2: (a) The dependence of the curvature-induced energy gap, -Bgap, on the chiral 
angle 9 and tube diameter dt- The surface is a plot of Eq. ^ which reproduces well 
the calculated results, (inset) Due to the nanotube curvature, the cutting line which was 
exactly crossing the Dirac point in the absence of curvature can miss the Dirac point 
when curvature is included. This curvature gives rise to an energy gap iJgap in "metallic" 
SWNTs. (b) The dt dependence of i?gap is given as a one-dimensional projection of 
(a) onto the dt axis. The points on the dashed, dot-dashed, and dotted curves satisfy 
n — m = 3, 6, 9, respectively. We plot the energy of the RBM, ^iwrbm of Eq. ([5]), as a solid 
curve for comparison. ^ 



gap 



assume that Ci = 223.5 [cm ^] and c?. = 12.5 [cm ^] which are e xperimentally 
derive d parameters as obtained by IStrano et al. (l2003h and iBachilo et al.l 
f l2002h . Using Eqs. © and (gD for zigzag SWNTs {6 = 0°), we see that 
is smaller than hw^BU when cit > 2[nm] (see Fig. Et^b)). 

The presence (absence) of a curvature-induced mini-energy gap in "metal- 
l ic" zi gzag (armchair) SWNTs was confirmed experimentally by lOuyang et al. 
( 1200 ll ). The chirality was measured experimentally for (9,0), (12,0), and 
(15, 0) zigzag SWNTs by these authors. The observed energy gap can be 
fitted by 4Ao / df which has the same dt dependence in Eq. (jl]) . Note that the 
coefficient is given by Aq = 37oa^j,/16, and 4Aq ^ 40[meV-nm^] is smaller 
than the value of c = 60[meV-nm^] in Eq. dl]). El This discrepancy may be at- 
tributed to ( 1) uniaxial and tor si onal strain which is un i ntentionally applied 
toaSWNT |Yang et al.l f ll999h : lYang and fl2000f ): iKleiner and Eggert 
(2001)], or (2) renormalization of the value of c due to the el-ph interaction, 
or (3) a SWNT-substrate interaction effect. (1,2) are intrinsic to SWNTs, 
while (3) is extrinsic. Since there are various factors which can affect the 
energy gap, it is not easy to predict the precise value of the energy gap, al- 
though the curvature-induced gap has been examined within the framework 
of first principles calcula tions including the effect of structure optimization 



Miyake and Saitd (120051 )]. It is noted that the chirality dependence of cos 3^ 
in Eq. @ has not been tested experimentally so far, except for 6 = (zigzag 
SWNTs) and 9 = 30° (armchair SWNTs). Study of a chiral SWNT is left 
for future experiments. 



3. Effect of Curvature on the Phonon Energy 



In this section we formulate the self-energy of a phonon mode (Sec. 13. ip . 
and explain qualitatively the effect of the curvature on the self-energy (Sec. l3.2( ). 



Putting 7o = 2.60[eV] and Occ — 0.142[nm] into the definition of Aq, we get tlie result 
AAq « 40[meV-nm2]. 

We expect tliat the curvature-induced gap follows (see Sec. 14.21 for the derivation) 

^^gap " I ^ ~ °} "^^^ 

when an uniaxial strain is applied to SWNTs. The value of a depends on the model used, 
but it is probably not dependent on dt- Considering the fact that the observed energy gap 
scales as d^^, the effect of strain is not so relevant. 



8 



The relationship between our formulation and that of others is referred to in 
Sec.lOl 



3.1. Phonon Self-Energy 

A renormalized phonon energy is written as a sum of the unrenormalized 
energy, fiw, and the real part of the self-energy, n(a;, E-p). The imaginary part 
of n(ci;, E-p) gives the spectrum width. Throughout this paper, we assume 
a constant value for for each phonon mode. The self-energy is given by 
time- dependent second-order perturbation theory as 

^(-'^-)-% ( fe._g'+.r/2 - iJT^TlrT^) >< < A " /=> • (^) 

where the pre-factor 2 comes from spin degeneracy, /h,e = (1 + exp(/3(£''^'^ — 
Ep))~^ is the Fermi distribution function, E^. {E^ is the energy of an electron 
(a hole) with momentum k, and E^ = E^^ — E^ (> 0) is the energy of 
an electron-hole pair. Vk is the el-ph matrix element that a phonon with 
momentum q = changes into an electron-hole pair [see the left diagram of 
Fig. [3](a)] which will be derived in Sec. HI Note that the momentum of an 
electron k is the same as that of a hole due to momentum conservation, and 
therefore pair creation involves a vertical transition. In Eq. ([7]), the energy 
shift is given by the real part of the self-energy, Re[n(c<;, i^p)], and the decay 
width r is determined self-consistently by r/2 = — Im [n(a;, ii^p)]- E The 
decay width relates to the average life-time r via r = h/T. It is noted that 
we use T = 300K although the self-energy is also a function of temperature 
= ksT where k-Q is Boltzmann's constant]. 

3.2. Phonon Softening and Hardening 

By defining the denominators of Eq. ([7]) as h±{E'^^) = ±l/{huj =]= E"^^ + 
ir/2), Eq. d?]) may be rewritten as 

U{u, Ep) = 2j2 W ) + h-{Et)] X (/h - fo) . (8) 

k 



^The self-consistent calculation begins by putting r/2 = 70 into the right-hand side of 
Eq. ([7]). By summing the right-hand side, we have a new r/2 via r/2 = — Im [n(w, iJp)] 
and we then put the new r/2 into the right-hand side again, iteratively. This calculation 
is repeated until n(aj, Ep) is converged. 
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When we assume that |Vkp does not depend on k, the E"^^ dependence of 
Re[U{u,E''^)] is determined by those of Re[h+{E^^)] and Re[/i_(E^^)]. It 
should be noted that Relh^^E"^^)] (sohd curve in Fig. Mjo)) has a positive 
(negative) value when E'^^ < hw {E'^^ > hw), and the lower (higher) energy 
electron-hole pair makes a positive (negative) contribution to Re[n(a;, ii^p)]- 
Therefore, the sign of the contribution to Re[n(ci;, Ep)], i.e., frequency hard- 
ening or softening, depends on its electron-hole virtual state energy, E'^^. 
In contrast, Re[h-{E'^^)] (dashed curve in Fig. [3](b)) always has a negative 
value, that is, it only contributes to a phonon softening. Note however that 
the contribution of Re[h-{E'^^)] is small compared with Re[/i+(£''^^)] since 
—l/hLO< Re[h_{E'^^)] < 0. Physically speaking, the h_{E'^^) term repre- 
sents an intermediate state including two phonons and electron-hole pairs 
(see the right hand diagram in Fig. [3](a)), while the h^^E'^^) term represents 
the intermediate state that includes only electron-hole pairs. |^ Even though 
the contribution of h-{E'^^) is relatively small, h-{E^^) is important to get 
a symmetric response of 11(0;, E-p) relative to the Fermi energy. In fact, due 
to the h_{E^^) term, the electron-hole pair at the Dirac point {E'^^ = 0) 
can not contribute to the self-energy, since Re[h^{E'^^) + h_{E'^^)] = when 
E""^ = 0. For high energy electron-hole pairs, the h±{E^^) terms contribute 
equally since Re[h+{E^^)] ^ Re[h_{E^^)] ^ -1/E^^. 

The curvature-induced energy gap, -Egap, affects the frequency shift since 
an electron- hole pair creation event is possible only when E^^ > -Egap- When 
< E'gap < hjj, the contribution to frequency hardening in Eq. ([7]) is sup- 
pressed. When ii^gap > hw, not only are all the positive contributions to the 
self-energy suppressed, but some negative contributions are also suppressed. 
Further, lm{h+{E^^)) is nonzero only when E is very close to hu, which 
shows that a phonon can resonantly decay into an electron-hole pair with the 
same energy. Thus, when Eg^^p > hw, we have F ~ because no resonant 
electron- hole pair excitation is allowed near E = fiw. It is therefore impor- 
tant to compare the values of -Egap and fkv for each (n, m) SWNT. For the 
LO and TO modes, hw is about 0.2 [eV] and therefore we get i?gap < fioj (see 
Fig. [2]) for most of the SWNTs except for a SWNT with a small diameter. 
Thus, those LO and TO modes can resonantly decay into an electron-hole 



^In fact, we have h.(E''^) = h+{2{nw + iT/2) + E""^). 

^For very small diameter SWNTs, the energy gap disappears because of the lowering 
of the interlayer energy bonds. 
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Figure 3: (a) In time-dependent second-order perturbation theory, we consider an inter- 
mediate state including only electron- hole pairs (the case oi t < i'), and an intermediate 
state including two phonons and electron- hole pairs (the case of t > t'). The former pro- 
cess corresponds to h+{E'^^), while the latter one corresponds to h^{E°^). (b) The energy 
correction to the phonon energy by an intermediate electron-hole pair, especially the sign 
of Jie{h+{E'^^)) (solid curve), that corresponds to frequency hardening or softening, de- 
pends on the energy of the intermediate state £'°'^. The contribution to n(w, E-p) of a low 
energy electron-hole pair satisfying < E^^ < E^^p is forbidden. Im(/i_|_(i?'^'')) (dashed 
curve) is nonzero only when £'°'^ is very close to fiw, which shows that a phonon mode can 
resonantly decay into an electron-hole pair with the same energy. 
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pair. The RBM mode in some SWNTs (for example, a (12, 0) zigzag SWNT) 
can not resonantly decay into an electron-hole p air, which re s ults in a long 



life-time for the RBM in that particular SWNT |Sasaki et all fl2008af )]. 

At T = 0, the Fermi distribution factor, namely /h — fe in Eq. (|7]), 
plays a very similar role as the curvature- induced gap, -Egap- In fact, all 
the excitations of electron-hole pairs with E*^^ < 2\Ep\ are forbidden due to 
the Pauli exclusion principle. A difference between the energy gap and the 
Fermi energy arises at a finite temperature. Some electron-hole pairs with 
-E'gap < < 2\Ef\ can contribute to the self-energy, while states E"^^ < E'gap 
do not exist even at a finite temperature. It should be noted that Vk in Eq. ([7]) 
depends on the value of E'gap since the position of the cutting line depends 
on i?gap, while Vk does not change by changing E^. This is also a crucial 
difference between the roles of -Egap and E-p in the self-energy. 

3.3. Other Formulas 

Here, we refer to the relationship between our formula and other formu- 
las. First, replacing r/2 in Eq. ([7]) with a positive infinitesimal 0+ gives the 
standard formula for the Fermi Golden rule. In this case, using l/(x-|-iO+) = 
P(l/x) — i7TS{x) with P denoting the principle value of integration and S{x) 
the Dirac delta-function, r/2 can be calculated directly, i.e., without using 
the self-consistent way, by performing the summation (or integral) of the 
right-hand side of Eq. ([7]). We calculate r/2 self-consistently by taking care 
of a finite energy level spacing originating from a finite length of a nan- 
otube where E^ now takes a discrete value, and is not a continuous variable. 
Roughly speaking, the broadening is suppressed when the energy level spac- 
ing, AE = 27ThvF/L, exceeds F. For example, the critical length where the 
broadening becomes negligible for a (10, 10) SWNT is about 700nm. 

Second, the summation index Ylk EQ- O is not restricted to only 
inter-band {E^^ ^ 0) processes but includes also intra-band (E"^"^ = 0) 
processes. Then, the self-energy can be decomposed into two parts, as 
n(w,EF) = Vl'-^^^^uj.Ey) + ff°*''^(w, Ef) where n'°*^'^(a;, Ep) includes only 
inter-band processes satisfying ^ 0. In the adiabatic limit, i.e., when 
w = and F = in Eq. ([7]), it is straightforward to get the following rela- 



®It may be appropriate to denote an intra-band process by an = or E = 
process. 
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tions, for a single Dirac cone at T = 0: 



U'''''''{0,Ef) 



k 



kl^X 



f(E^),(EF>0) 

f(i?^), (Ep<0) 



--\E.\ 
2' ' 



VIP 



a 



(9) 



k 



-2 (^^-I^fI 



where a = S\v\'^/TT{hvpy and Ec is some cut-off energy. Here, we have 
assumed that = v cos e(k). E Note that ff^*''^(0, Ep) does not vanish 
because (/h — fe)/E^ 7^ in this hmit, while in the non-adiabatic case, 
n™*™(a;, iJp) vanishes since (/h — fe)/hu = 0. It is only the inter-band 
process that contributes to the self-energy in the non-adiabatic case. In 
the non-adiabatic limit at T = 0, it is a straightforward calculation to get 



a 



Re[U{uj,EF)] = -- 



E, 



EJ -In 



\E-p 


2 




2 



, (E, >M, (10) 



where E^ 



2hv-pk, V / {2'kY j^" kdk Jq"" dQ, and / x/{x + a)dx = 

X — aln\x + a\ have been used in Eq. ([7]) to get the right-hand side. The 
Fermi energy dependence is given by the last two terms for the case of a 
massless Dirac cone spectrum. The first term is linear with respect to Ep 
and the second term produces a singularity at \Ep\ = hw/2. This singularity 
is useful in identifying the actual Fermi energy of a graphene sample. 

It is also interesting to consider the case of a massive Dirac cone spectrum, 
E = vn? + {hypky. In the non-adiabatic limit at T = 0, we get 

{huf - (2m)2 



a 



Re[U{u,Ep)] = -- 



Er — Ep — 



In 



I^fI 


tlLU 

2 


I^fI 


2 



where Ec ^ hio and = v{hvpk/E) cosQ(k) are assumed. Equation ( ITTjl 
is for I^'fI > m. For \Ep\ < m, the self-energy shift is given by replacing 



10 



8(k) is the angle between the vector k and the fc^j-axis [see Eq. (fT4] 



Lazzeri and Mauril ( 2006 ) showed that 11(0, iJp) does not depend on Ep in the adi- 
abatic hmit due to the canceUation between n*"*'''^(0, £;f) and n''^*'=''(0, E'p). This shows 
that the adiabatic approximation is not appropriate for discussing the Ep dependence of 
the self-energy. 
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I E-p I with m in Eq. fllip . The logarithmic singularity for the last term disap- 
pears when m = hLo/2, and its overall sign is interchanged when 2m > hu. 
Broadening is possible only when 2m < hw, which may be useful in knowing 
whether the graphene sample has an energy gap or not. 



4. The Electron-Phonon Interaction 

In this section we provide a framework to obtain the el-ph (electron- 
phonon) interaction in the effective-mass theory, and show how to calculate 
the el-ph matrix elements. The main results are Eqs. (!22|) and ( 147|) . Those 
who are not interested in the details of the derivation can skip this section. 

4.1. Unperturbed Hamiltonian 

The unperturbed Hamiltonian in the effective-mass model for vr-electrons 
near the K point of a graphene sheet is given by 

= VYO- ■ p, (12) 

where vy is the Fermi velocity, p = —iKV is the momentum operator, and 
cr = {(7x,(7y) is the Pauli matrix. The x, y, and z coordinate system is 
taken as shown in Fig. |H is a 2 x 2 matrix which operates on the two 
component wavefunction: 

■*^(^) - (M) ■ 

where V'a(^) V'b('') ^'^^ wavefunctions of vr-electrons for the sublat- 
tices A and B, respectively, around the K point. The energy eigenvalue of 
Eq. f ll2p is given by ±fF|p| and the energy dispersion relation shows a linear 
dependence at the Fermi point, which forms what is known as the Dirac cone. 

The energy eigenstate with wave vector k in the conduction energy band 
is written by a plane wave e^^'^ with the Bloch function ipf^^ as ipfi^{r) = 



We use the Pauli matrices of the form of ax ~ o) ' ~ i^i 0^^ ' ^^'^ '^^ ~ 
Q ) ' "^^^ 2x2 identity matrix ctq is given by cto = f g 'j' 
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j^gik r^K^ where iV is a normalization constant satisfying N'^S = 1, 5" is the 
area (volume) of the system, and 

<^ ^ (ei(k)) • (14) 

Here k is measured from the K point, and ©(k) is defined by an angle of 
k = {kx,ky) measured from the kx-3.xis as {kj.,ky) = |k|(cos0(k), sin 9(k)). 
The eigen value of this state is = +?;f|p|- The energy eigenstate with the 
energy eigen value E = —vf\p\ in the valence energy band is written by 

<k(r) = ^ [_J^ei^)) ■ (15) 

The energy eigenstate for the valence band, ^^k(r) is given by azip^i^ir). This 
results from the particle-hole symmetry of the Hamiltonian: (TzI-L^cFz — —'H^- 




Figure 4: A hexagonal unit cell of graphene consists of A (closed circle) and B (open 
circle) sublattices. ai and a2 are lattice vectors. Ra (a = 1,2,3) are vectors pointing to 
the nearest-neighbor B sites from an A site (Ri = flccGj,, R2 = — (\/3/2)(iccea; — (l/2)accey, 
and R3 = (V3/2)accea; — (l/2)accey) where (ey) is the dimensionless unit vector for 
the X-axis (y-axis). Local modulations of the hopping integral are defined by (r^) 
(a = 1, 2, 3) where is the position of an A- atom. 

The unperturbed Hamiltonian near the K' point is given by 

nf = v^a' ■ p, (16) 
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where cr' = (— cr^, (Jy). The dynamics of vr-electrons near the K' point relates 
to the electrons n e ar th e K point by time-reversal symmetry, ip^ — )■ {ip^ )* 
Sasaki and Saitol fl2008h ]. Because lattice vibrations do not break time- 



reversal symmetry, we mainly consider electrons near the K point in this 
paper. 

4-2. Deformation- Induced Gauge Field 

Lattice deformation modifies the nearest-neighbor hopping integral lo- 
cally as —70 — )■ —70 + 57o(ri) (a = 1,2,3) (see Fig. H]). The corresponding 
perturbation of the lattice deformation is given by 

^1 = E E ^^oir.) [(cf+J^cf + (cf )tcf^,] , (17) 

ieA a=l,2,3 

where cf is the annihilation operator of a vr electron of an A-atom at position 
Tj, and {cf^aV is a creation operator at position rj+a (= + Ra) of a B-atom 
where Rq (a = 1,2,3) are vectors pointing to the three nearest-neighbor B 
sites from an A site. 

The perturbation of Eq. ([T7j) gives rise to scattering within a region 
near the K point (intravalley scattering) whose interaction is given by a 
deformation-induced gauge field A(r) = {Ax{r), Ay{r)) in Eq. f|T2l) as 

n^ + nf = vpa-[p + A\r)], (18) 



wher e A'^fr) is defined from^7n,a (r) (a = 1,2,3) as [Sasaki et al.l (12005 



20061 ): IKatsnelson and GeimI (1200 



VpAlir) = 570,1 (r) - i [570,2 (r) + 57o,3(r)] 
z;FA^(r) = ^[57o,2(r)-57o,3(r)]. 



(19) 



When 570,2 = 570,3 = 0, then A%r) = (A^(r),0) and Ai(r) ■ Ri = 0. 
Similarly, when 570,1 = 570,3 = 0, we have A'i(r) ■ R2 = 0. Generally, 
the direction of A*^(r) is pointing perpendicular to the bond whose hopping 
integral is changed from 70. For the K' point, we obtain 

nf +nf' =iW -Ip-A^r)]. (20) 

Even though the A^(r) appears as a gauge field, it does not break time- 
reversal symmetry because the sign in front of A'^(r) is opposite to each 
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other for the K and K' points. This is in contrast with the fact that A(r) 
(vector potential) violates time-reversal symmetry because the sign in front 
of A(r) is the same for the K and K' points since p — p — eA(r) in the 
presence of a magnetic field. 

The gauge field description for the lattice deformation (Eq. f|T9|) ) is useful 
to show the appearance of the curvature- induced mini-energy gap in metallic 
carbon nanotubes. For a zigzag nanotube, we have 570,1 = and ^70,2 = 
(570,3 7^ from the rotational symmetry around the tube axis (see Fig. H]). 
Then, Eq. (fT9l) shows that for 7^ and A'^ = 0, the cutting line of = 
for the metallic zigzag nanotube is shifted by a finite constant value of A"^ 
because of the Aharanov-Bohm effect for the lattice distortion-induced gauge 
field A^. For an armchair nanotube, we have ^70,1 7^ and ^70,2 = ^10,3- 
Then, Eq. f|T9|) shows that for A'^ ^ and = 0, the cutting line of ky = 
for the armchair nanotube is not shifted by a vanishing A^. This explains 
the presence (absence) of the curvat ure-induced mini - energ y gap in metallic 



zigzag (armchair) carbon nanotubes Kane and Meld (119971 )] 



The gauge field description is also useful to discuss the effect of an uniaxial 
strain on the gap. Let us consider applying a strain along the axis of a zigzag 
SWNT. Then, due to the symmetry, we have ^70,1 = a and ^70,2 = ^70,3 = 
a/2 where a is a constant. Putting these perturbations into Eq. (fT9|l we see 
that vpA"^ = a/2, which means that the curvature- induced gap in a zigzag 
nanotube can change a little by the strain along the axis. For an armchair 
SWNT, instead, we have ^70,1 = and 570,2 = ^70,3 = b, which results in 
vpA"^ = 0. This shows that the absence of the gap in armchair SWNT is 
robust against a strain applied along the nanotube axis. 

4-3. Deformation- Induced Gauge Fields for LO and TO modes 

Here, we derive A'^(r) for the LO and TO modes. Let u(r) be the relative 
displacement vector of a B site from an A site (u(r) = UB(r) — UA(r)) and 
let g be the el-ph coupling constant, then 57o,a(r) for the LO and TO modes 
is given by 

<57o,a(r) = |u(r) • (21) 

where Ha denotes the nearest-neighbor vectors (Fig. [H and Fi g . [Ma) ) and 
g = 6.4 eV/A is the off-site el-ph matrix element Porezag et al. ( 19951 )]. We 
rewrite Eq. ( !T9|) as 

Vp{Al{r),A^(r)) = g{uy{r),-u,{r)), (22) 
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where Ui{r) = u(r)-ej, {i = x,y), and Ri — (R2+R3)/2 = iey and -\/3/2(R2 — 
R3) = —ie^ {i = 3acc/2) have been used (see the caption of Fig. Hj). Then, 
the el-ph interaction for an in-plane lattice distortion u( r ) can be rewritten 
as the vector product of cr and u(r) jlshikawa and Ando ( 2006 )] as 



Kg = vfo- ■ Ai(r) = g{(T x u(r)) ■ e. 



(23) 




Figure 5: (a) The hexagonal lattice deformed by a phonon displacement vector. Let 
UA(r) (uB(r)) is a displacement vector of an A (B) site, then the modulation given by 
optical phonon modes is (57o,a = (ff/^)u(r) • Rq where u(r) {— UA(r) — UB(r)) is a relative 
displacement vector of a B site relative to the nearest A site, (b) The cross section of 
a nanotube. The displacement vector for the RBM, s(x), is decomposed in terms of the 
normal Sz{x) and tangential Sx{x) components. The derivative of the normal unit vector 
Bz (x) with respect to x gives a component along e^; (x) , which modifies the net displacement 
along the x direction. 

The gauge field description for the el-ph interaction of the LO and TO 
modes (Eq. ( l22i) ) is useful to show the absence of the el-ph interaction for the 
TO mode with a finite wavevector, as shown below. The TO phonon mode 
with q 7^ does not change the area of the hexagonal lattice but instead 
gives rise to a shear deformation. Thus, the TO mode (uTo(r)) satisfies 

V-UTo(r) = 0, VxuTo(r)^0. (24) 

Using Eqs. (122!) and (l24l) . we see that the TO mode does not yield a deformation- 
induced magnetic field, 

B^(r) = V X A^(r), (25) 
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but the divergence of A''(r) instead does not vanish because 
5Kr) = -^V-UTo(r)=0, 

(26) 

Vf 

Thus, we can define a scalar function (f{r) which satisfies A''(r) = Vip{r). 
Since we can set A^( r) = in Eg . (fTSll b v selecting the gauge as ip^^r) — >■ 



exp{-iip{r)/h)ip^{r) [Sasaki et all (120051 )] and thus the A^(r) in Eq. 



disappears for the TO mode with q 7^ 0. This explains why the TO mode 
with q 7^ is completely decoupled from the electrons, and that only the 
TO mode with q = couples with electrons. This conclusion is valid even 
when the graphene sheet has a static surface deformation. In this sense, the 
TO phonon mode at the F-point is anomalous since the el-ph interaction 
for the TO mode can not be eliminated by a phase of the wavefunction. In 
contrast, the LO phonon mode with q 7^ changes the area of the hexagonal 
lattice while it does not give rise to a shear deformation. Thus, the LO mode 
(uLo(r)) satisfies 

V-ULo(r)^0, VxuLo(r) = 0. (27) 

Using Eqs. fl22l) and fl27|) . we see that the LO mode gives rise to a deformation- 
induced magnetic field since 

5,(r) ^0, V ■ A(r) = 0. (28) 

Since a magnetic field changes the energy band structure of electrons, the 
LO mode can couple strongly to the electrons even for q 7^ 0. 

4-4- Deformation- Induced Gauge Field for the RBM 

Next, we derive the deformation-induced gauge field A^(r) for the RBM. 
When the RBM displacement vector of a carbon atom at r is s(r) = (sa;(r), Sy{v), ^^(r)), 
the perturbation to the nearest-neighbor hopping integral is given by 

<57o,a(r) = |r„ ■ {s(r + R,) - s(r)}, (29) 

By expanding s(r + R^) in a Taylor's series around the displacement s(r) as 
s(r + Ra) = s(r) + (R^ ■ V)s(r) + ■ ■ ■ , we approximate Eq. (12^ as 

57o,a(r)~fRa-{(Ra-V)s(r)}. (30) 
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Putting Ri 



-(v^/2)c 



(l/2)accej^, and R3 = {\/?,/2)accGx- 



(1/2) 

e^, into the right-hand side of Eq. (|30|), we obtain the corresponding 
deformation-induced gauge field of Eq. f lT9|) as 



v^Alir) 



gac 



2 

ga-o 



ds:. 



dx 
dsJr) 



dy 



+ 



^ dsy jr) 
dy 

dx 



(31) 



Further, the displacements of carbon atoms give an on-site deformation po- 
tential in which the diagonal Hamiltoni a n matrix elements are modi fied by 
the el-ph interaction [Jiang et al.l ( 120051 ): ISaito and Kamimural (Il983[ )] 



fi'onO'O 



dsjj) 
dx 



+ 



dsy jr) 
dy 



(32) 



Here, dsx{'r)/dx + ds JT)ldy {= V • Sfr) ) repre sents the change of the area 
of a graphene sheet Suzuura and Andq ( 2002 )]. According to the density 
functional calculation by Porezag et al.~fll995l ). we adopt the on-site coupling 
constant = 17[eV]. 

Since Eqs. (13T!) and (!32|) are proportional to the derivatives of Sx{r) and 
Sy{r), that is, they are proportional to q, where q is the phonon wave vec- 
tor. Then, the el-ph matrix element for the in-plane longitudinal/transverse 
acoustic (LA/TA) phonon modes vanishes at the P point Namely, A''(r) = 
and Hon = in the limit of q = 0. Among the TA phonon modes, there 
is an out-of-plane TA (oTA) phonon mode. The oTA mode thus shifts car- 
bon atoms on the flat 2D graphene sheet in the z-direction [see Fig. H] and 
Fig. El^b)]. The oTA mode of graphene corresponds to the RE M of a nan- 
otube even though the REM is not an acoustic phonon mode Saito et al 
(119981 )]. In the following, we will show that the el-ph interaction for the 
REM is enhanced due to the curvature of the nanotube as compared with 
the oTA mode of graphene since the REM is a bond-stretching mode due to 
the cylindrical structure of SWNTs. 

The displacements of the REM modify the radius of a nanotube as r ^ 
r + Sz{r) (see Fig. [5](b)). A change of the radius gives rise to two effects to the 
electronic state. One effect is a shift of the quantized transverse wave vector 
around the tube axis. The distance between two wave vectors around the 
tube axis depends on the inverse of the radius due to the periodic boundary 
condition, and a change of the radius results in a shift of the wavevector. 
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The other effect is that the RBM can change the area on the surface of 
the nanotube even at the F point. This results in an enhancement of the 
on-site el-ph interaction. These two effects are relevant to the fact that the 
normal vector on the surface of a nanotube is pointing in a different direction 
depending on the atom position. To show this, we take a (zigzag) nanotube 
as shown in Fig. \5Kh)- Let us denote the displacement vectors of two carbon 
atoms at x and x + dx as s{x) and s{x + dx), then an effective length for the 
displacement along the between the nearest two atoms is given by 

Dx = ex{x + dx) ■ [s{x + dx) — s(x)] . (33) 

By decomposing s{x) in terms of a normal and a tangential unit vector as 
s{x) = Sz{x)ez{x) + Sx{x)ex{x) (see Fig. ID^b)), we see that Eq. (133|1 becomes 

Dx = Sx{x + dx) + Sz{x + dx)ex{x + dx) ■ e^^x + dx) 
- Sx{x)ex{x + dx) ■ ex{x) - Sz{x)ex{x + dx) ■ e^^x) 
, ( dsx(x) Sz(x)^ 



dx 

where we have used the following equations: 

ez{x + dx) = ez{x) H e^.(x) + 

r 

ex(x + dx) = ex(x) H e^(x) + 

r 



(35) 



Equation (1341) shows that the net displacement along the x axis is modified by 
the curvature of the nanotube as dxSx{r) — )■ dxSx{r) + Sz{r)/r. The correction 
is negligible for a graphene sheet (r — )■ oo), but appears as an enhancement 
factor to the el-ph interaction in SWNTs. 

The el-ph interaction for the RBM is included by replacing dxSx{r) with 
dxSxir) + Sz(v)/r in Eqs. fl3T|) and fl32l) . In Eq. fl3T|) . we have an additional 
deformation-induced gauge field, 

v^Al{v) = -^ffM, (36) 

for the RBM mode which gives rise to a shift of the wavevector around the 
tube axis even at q = 0. In Eq. (152]) . it is shown that the RBM produces an 
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additional on-site deformation potential of Qon'^oiszi'^)/''^)- Finally, we obtain 
the el-ph interaction for the F point (q = 0: s(r) is a constant) RBM, as 



HRBM = -^-a.. + ^7o„-^o = ^^_^ ). (37) 

This representation is for zigzag SWNTs. For a general (n, m) SWNT with 
a chiral angle 9, the el-ph interaction for the RBM becomes 

^rbmI^*) = I _ga,Tp-m ^„ • (38) 

Ut \ 2 / 



See ISasaki et al.l ( l2008al ) for more details. 



5. Kohn Anomaly Effect 

Here we consider the el-ph matrix element as a function of the electron 
wavevector k for the LO and TO phonon modes and the RBM with q = 
(i.e., F-point). The displacement vector with q = is expressed by a position 
independent u = {ux,Uy), by which an electron-hole pair is excited. The el- 
ph interaction with q = is relevant to phonon-softening phenomena for all 
three kinds of modes. 

5.1. Matrix Element for Electron-hole Pair Creation 

Let us first consider the case of a zigzag SWNT. In Fig. HI we denote y 
(x) as a coordinate along (around) the axis of a zigzag SWNT, and Uy (ux) 
are assigned to the LO (TO) phonon mode. Thus, from Eq. fl22l) . we have 

vfAIo = g{uy,0), ^^^^ 
vpA'^Q = g{0, -Ux). 

The direction of the gauge field A^(r) is perpendicular to the phonon eigen- 
vector u and the LO mode shifts the wavevector around the tube axis, which 
explains how the LO mode may induce a dynamical energy band-gap in 



^^In case of the F point phonon, the definition of the LO and TO is not unique. It seems 
standard that the LO is taken as the mode parallel with respect to the tube axis and the 
TO mode is the one perpendicular to the tube axis. 
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metallic nanotubes |Dubay et al.l (|2002[ )]. Putting Eq. ( 15^ into Eq. (12^ . we 
get 



nil = vyKo ■ ^ = 9Uya^, .^Q. 

The el-ph matrix element for the electron-hole pair generation is given 
from Eqs. dS]), and (00]), by (A = LO,TO) 

{eHk)\Ht\oo,) ^ J (41) 

By calculating Eq. (HT]) for the LO mode with {ux,Uy) = (0,m) and for the 
TO mode with {ux,Uy) = {u,0), we get 

(eh(k)l'HLol^Lo) = -igusmOik), ^^^^ 
(eh(k)|7{^^^|(X'To) = — igu cos Q(k), 

where 0(k) is defined by an angle of k = {kx, ky) measured from the k^ axis. 

Next, we consider the case of an armchair SWNT. In Fig. HI x (y) is the 
coordinate along (around) the axis and Ux {uy) is assigned to the LO (TO) 
phonon mode. Then, for an armchair SWNT, we get 

{ehiX)\HZ^\ui^o) = -^gusm^{k), 

{eh{k)\H^'^\ujTo) = -igu cos e{k). ^ ^ 

Note that 6(k) for the armchair nanotube is given by rotating G(k) for the 
zigzag nanotube by 7r/2 (^^(k) = 6(k) + tt/2). It is useful to define the ki 
(^2) axis pointing in the direction of a general SWNT circumferential (axis) 
direction (see Fig. [6]), and d(k) as the angle for the polar coordinate. Then, 

(eh(k)lHLol^^Lo) = -i5'^^sin6'(k), ^^^^ 
(eh(k)|'HTo|i^To) = — igu cos 6(k). 

is valid regardless of the tube chirality if the phonon eigenvector of the LO 
(TO) phonon mode is in the direction along (around) the tube axis. This is 
because p and u(r) [a nd A^(r)] are transfo rmed in the same way when we 
change the chiral angle [Sasaki et al. fl2008ah ]. As a result, there would be no 



chiral angle dependence for the el-ph matrix elements in Eq. (jH]). Note also 
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Figure 6: Dependence of LO and TO phonons on the cutting line, (a) A cutting line near 
the K-point. The fci (fc2) axis is selected as the nanotube circumferential (axis) direction. 
The amplitude for an electron-hole pair creation depends strongly on the relative position 
of the cutting line from the K-point. (b) If the cutting line crosses the K-point, then the 
angle 9{k) (= arctan(fc2/fci)) takes tt/2 (— 7r/2) values for fc2 > (fc2 < 0). In this case, 
the LO mode strongly couples to an electron-hole pair, while the TO mode is decoupled 
from the electron- hole pair according to Eq. (|43p . 
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that Eq. f Hl|) shows that {eh.(k)\'Hx\oJx) depends only on 9(k) but not on |k|, 
which means that the dependence of this matrix element on E'^^ (= 2/iVF|k|) 
is negligible [see Fig. [6](b)]. 

Where does the ^(k) dependence in Eq. (jH]) then come form? The expec- 
tation value of ax, (Ty, and az with respect to ^'^^^(r) defines the pseudospin. 
Using Eq. f[T^ with 6(k) — t- ^^(k), we have the expectation values for the 
Pauh matrices (a^.) = {tpf^kWxH'f^i,) = cos6'(k), (a^.) = (^A^kl^yl^^^k) = 
sin6'(k), and (ax) = (i'^-^Wzli'^k) = 0. Then the direction of the pseudospin 

ofV'c^W, 

{{ax), {ay), {a,)) = (cos ^(k), sin ^(k), 0), (45) 

is within {ki, ^2) plane and parallel to k (see Fig. [7]). Due to the particle- 
hole symmetry, ip^y_{r) = az'ip^]^{r), the el-ph matrix element for the electron- 
hole pair creation process can be related to the pseudospin. For example, we 
see that 

(eh(k)|Kg|a;Lo) = 9Uy{^lM^^,^) 
= -igu{ay) 

= — igu sin 9{k). (46) 

The electron-hole pair creation for the LO mode is relevant to the pseudospin 
component which is parallel to the tube axis, {ay), while that for the TO 
mode is relevant to {a^). 

For RBM, from Eq. (!38l) . the matrix element for an electron-hole pair 
creation is chirality dependent as 

(eh(k)|HRBM(^)|c^RBM) = ^ga^c^ sm{e{k) + 3d). (47) 

Thus, the frequency shift of the RBM can have a chiral angle dependence. In 
particular, armchair SWNTs {6 = 30°) exhibit neither a frequency shift nor 
a broadening, regardless of their diameters because the el-ph matrix element 
becomes 

(eh(k)|?^^^^Ml^RBM) = tga,,^ cose{k), (48) 
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For the pseudospin of the electrons near the K' point, see ISasaki et ahl (|201 
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Dirac cone ^2 

Pseudospin 




Figure 7: The el-ph matrix element for the electron-hole pair creation process relates to 
the pseudospin of the electronic state. 



which is zero for a cutting hne for a metaUic band: ^(k) = ±7r/2. This 
^(k) dependence of Eq. f HHj) is the same as that of the TO phonon mode of 
Eq. fH3|) . so that the absence of a frequency shift of the RBM in armchair 
SWNTs is similar to the abse nce of a frequency sh ift of the TO mode at the 
r point in armchair SWNTs Sasaki et al.l ( l2008bh ]. 



5.2. Phonon frequency shift 

Here we show the calculated results for the phonon frequency as a function 
of the Fermi energy. 

5.2.1. Armchair SWNTs 

First, we consider Eq. ( H4l) for a fc-point (k = (A;i, ^2)) on the cutting line 
of an armchair SWNT. Since the armchair SWNT is free from the curvature 
effect, the cutting line for its metallic energy band satisfies fci = and lies 
on the /c2 axis. Thus, we have ^(k) = 7r/2 (— 7r/2) for A;2 > {k2 < 0). Then, 
Eq. (jH]) tells us that only the LO mode can couple to an electron-hole pair 
and the TO mode does not couple to an electron-hole pair for the metallic 
energy band of an armchair SWNT. Similarly, Eq. (HHj) shows that the RBM 
of an armchair SWNT does not show any phonon softening. 

In Fig. [HI we show the phonon energy as a function of E-p for a (10, 10) 
armchair SWNT. Here we take 1620 cm~^ and 1590 cm~^ for huj of the LO 
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and TO modes, respectively. The energy bars denote T values. The self- 
energy is calculated for T = 300K and L = 10/im. It is shown that the TO 
mode does not exhibit any energy change, while the LO mode shows both an 
energy shift and a broadening. As we have mentioned, the minimum energy 
is realized at \Ep\ = fvjj/2 0.1 eV). There is a local maximum for the 
spectral peak at \Ep\ = 0. The broadening for the LO mode has a tail at 
room temperature for \E-p\ > hw/2. 




Ef (eV) 



Figure 8: (color online) The Ep dependence of the LO (red curve) and TO (black curve) 
phonon energy in the case of the (10, 10) armchair SWNT. The data are taken at room 
temperature. Only the energy of the LO mode is shifted, with the TO mode frequency 
being independent of Ep. The decay width (F) is plotted as an error-bar. The blue curve 
is given by the analytic result of Eq. . 



In evaluating the LO mode's self-energy according to Eq. ([7]), we have 
assumed that the cutoff energy is Ec = 0.5eV. The presence of a cutoff energy 
is reasonable since the matrix element act ually depends on the energy of the 
electron-hole pair [see ISasaki et al.l (|2009[ )]. An analytical expression for the 
Ec dependence of the self-energy is easy to obtain by using the effective-mass 
model, which can be derived from Eq. ([7]) at T = as 



Re [n(w, Ep)] 



L 



TT 



h 



In 



+ In 



\Ep\ + %^ 



(49) 
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The factor scales as \/dt because 

For comparison, we plot /k^; + Re [n(cj, E-p)] (Eq. (149|) ) as the blue curve in 
Fig. [HI 

5.2.2. Zigzag SWNTs 

Next we consider "metallic" zigzag SWNTs. When the curvature effect 
is taken into account, the cutting line does not lie on the K-point, but is 
shifted by ki from the ^2 axis. In this case, cos6'(k) = ki/{k\ + /c|)^/^ is 
nonzero for the lower energy intermediate electron-hole pair states. Thus, 
the TO mode can couple to the low energy electron-hole pair which makes 
a positive energy contribution to the phonon energy shift. The high energy 
electron-hole pair still decouples from the TO mode since cos0(k) — t- for 
\k2\ > \ki\. 

In Fig. [U^a), we show calculated results for the LO and TO modes as a 
function of for a (12, 0) zigzag SWNT. In the case of zigzag SWNTs, not 
only the LO mode but also the TO mode couples with electron-hole pairs. 
The spectrum peak position for the TO mode becomes harder (upshifted) for 
E-p = 0, since Re{h+{E'^^)) for E"^^ < tvjj contributes to a positive frequency 
shift. The hardening of the TO mode is a signature of the curvature- induced 
mini-energy gap. 

In Fig. El^b), we show the result for the RBM. The matrix element of 
Eq. fH7|) with 6* = 0° is proportional to sin6'(k). Thus, the high energy 
electron-hole pair can couple to the RBM and can contribute to the softening 
of the RBM. Although the magnitude of the shift is smaller than those for 
the LO mode , the softening for the RBM can be observed experimentally 



Farhat et al.l (120091)]. 



5.2.3. CUral SWNTs 

Finally, we examine "metallic" chiral SWNTs. The same discussion for 
the "metallic" zigzag SWNTs can be applied to "metallic" chiral SWNTs. 
However, there is a complication specific to chiral SWNTs that the phonon 



Here we use a harmonic oscillator model which gives u — y/h/2MJ^^ where Mc is 
the mass of a carbon atom. Using fiio — 0.2eV, we get y/Nu\A'^/h\ ss 2 x 10~^A~^. 
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Figure 9: (color online) (a) The Ep dependence of the LO (red curve) and TO (black 
curve) phonon energy in the case of the (12, 0) zigzag SWNT. The data are taken at room 
temperature. Not only the frequency of the LO mode but also that of the TO mode is 
shifted due to the curvature effect. Here we take 1640 cm~^ and 1590 cm~^ for hui of the 
LO and TO modes, respectively, (b) The Ep dependence of the REM energy hu! in the 
case of the (12,0) zigzag SWNT. 



eigenvector depends on the chiral angle. iReich et al.l (120011 ) reported that, for 
a chiral nanotube, the atoms vibrate along the direction of the carbon-carbon 
bonds and not along the axis or the circumference. The phonon eigenvector 
of a chiral nanotube may be written as 



Uto 

ulo 



COS (p sm ( 
— sin (h cos < 



(51) 



where Ui (^2) is in the direction around (along) a chiral tube axis, and (j) 
is the angle difference between the axis and the vibration. This modifies 
Eq. (03]) as 

(eh(k)l'HLolwLo) = -igusm{e{k) + 0), 
(eh(k)l'HTolc^To) = -^5'^^ cos(6'(k) + 0). 

The identification of in Eq. (I52p as a function of chirality would be useful 
to compare theoretical result s and experiments, which will be explored in the 
future [See iPark et al.l ( 120091 ) for example] . 



(52) 



5.2.4- Graphene 

In the case of 2D graphene, Eq. (j^3|) tells us that the F point TO and 
LO modes give the same energy shift because the integral over 6(k) gives 
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the same self-energy in Eq. ([7]) for both TO and LO modes. Th is explains 



why n o G-band splitting is observed in a single layer of graphene |Yan et al. 



(120071 )]. Even when we consider the TO and LO modes not exactly at the 
r point, we do not expect any splitting between the LO and TO phonon 
energies since the TO mode with q 7^ is completely decoupled from the 
electrons [see Eq. f l2B]) ]. Thus, for q 7^ 0, only the LO mode contributes to 
the G band intensity. It is interesting to note that the F point LO and TO 
modes may exhibit anomalous behavior near the edge of graphene because 
the wave function is not given by a plane wave but rather by a standing wave. 
The pseudospin for the standing wave is different from that for a plane wave. 
Moreover, the standing wave near the zigzag edge is different from that near 
the armchair e dge, which gives rise to a selection rule in their Raman spectra 



Sasaki et al.l ( 120091 . |2010[ )]. The standing wave behavior near the edges in 



graphene ribbons is beyond the scope of the present paper. 

6. Discussion and Summary 

We have seen that the curvature- induced gap is absent for armchair 
SWNTs, so that here the LO mode exhibits a strong Kohn anomaly ef- 
fect. Recently, however, it has been reported that even armchair SWNTs 



have an energy gap originating from a correlation effect [Deshpande et al. 



( I2OO9I )]. The correlation-induced gap observed is approximately 80 meV for 
armchair SWNT with dt = 2nm, and the gap increases with decreasing dt- 
Since the presence of a gap suppresses the contribution to the hardening, a 
local maximum (around E-p = 0) in the LO frequency vs. Ep plot [see Fig. |8] 
may disappear, and E-p = would become a global minimum if the correla- 
tion gap exceeds 200meV. Confirming this behavior in a Raman spectroscopy 
study may provide further evidence for the correlation-induced gap. 

Finally, we discuss our results in relation to the experimental results of the 



Kohn anomaly for the RBMs in metallic SWNTs |Farhat et al.l f l2009f )]. The 
frequency shifts observed are approximately 2cm~^, which is smaller than 
the theoretical value for a (12,0) zigzag SWNT, 8cm~^, shown in Fig. lU^b). 
This discrepancy may be attributed to the choice of the off-site el-ph matrix 
element, g, although we have determined this value from density functional 
calculation. Indeed, by decreasing the value of g from g = QA eV/A we find 
a better agreement to the experimental result, 3cm~^ when (7 = 4 eV/A. 
Another possibility is that we have used a harmonic oscillator model for 
obtaining the magnitude of displacement vector u. Actual value of u may 
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be smaller than our estimation, which also gives a better agreement to the 
experimental result. 

In summary, the el-ph interaction with respect to the LO, TO, and RBM 
for the Raman feature for SWNTs is derived in a unified way using the 
deformation-induced gauge field. Then, we have shown that the matrix ele- 
ment for electron-hole pair creation depends on the position of the cutting 
line. As a result, the TO mode in "metallic" SWNTs, except for armchair 
SWNTs, can couple to an electron-hole pair due to the curvature effect which 
shifts the cutting line away from the K point. In particular, only the low 
energy electron-hole pairs can couple to the TO mode and give rise to a 
hardening of the TO mode. The hardening of the TO mode is suppressed 
for large diameter SWNTs. This is reasonable since the TO mode as well as 
the LO mode exhibit a softening in the case of graphene samples. 
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